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Earth's magnetic field is essential for orientation in birds migration (T}|6l. The most promising explanation 
for this orientation is the photo- stimulated radical pair (RP) mechanism fT-TTl, conjectured to occur in cryp- 
tochrome photoreceptors 1 12 13 1. The radicals must have an intrinsic anisotropy in order to define a reference 
frame for the compass 1 14- 16|. This anisotropy, when introduced through hyperfine interactions, imposes im- 
mobility of the RP, and implies that entanglement between the unpaired electrons of the RP is preserved over 
long coherence times ifTTllTSl . We show that this kind of anisotropy, and consequently the entanglement in the 
model, are not necessary for the proper functioning of the compass. Classically correlated initial conditions for 
the RP, subjected to a fast decoherence process, are able to provide the anisotropy required. The environment in 
which the RP is immersed is then responsible for the reference frame of the compass, relaxing the immobility 
assumption. This fact significantly expands the range of applicability of the RP mechanism providing more 
elements for experimental search. 



The ability of birds to use the Earth's magnetic field to ori- 
entate themselves in the correct direction for migration (T] [2| 
has originated several experimental works devoted to the un- 
derstanding of the main features of the underlying mechanism 
(SHSI. One of the first proposals for modeling the appearance 
of the magnetic compass was that magneto-perception oper- 
ates by means of anisotropic magnetic field effects on the rate 
of production of yields of a photo stimulated radical pair re- 
action IMHEI . Although other models have been proposed, 
such as a magnetite-based magneto-perception [19, 20 1, the 
proposal of a radical pair mechanism (RPM) has been recently 
reinforced |7 1 and strong experimental evidences have been 
presented in its favor |8-10, 2ll|22l. The solid state RPM 
model can be summarized as follows |7|: a molecular precur- 
sor reacts to form a pair of radicals due to photochemically- 
driven electron transfer. Taking into account that both radi- 
cals are created in a single event, it is natural to assume that 
the electron spins are initially entangled (in the following we 
assume that the radical pair (RP) is created in a singlet s state, 
although working with a triplet t state is also possible). This 
singlet state evolves under the influence of a Hamiltonian con- 
taining an hyperfine interaction term between the nuclei and 
their electrons and a Zeeman interaction term between the un- 
paired electrons. Due to the anisotropy of the hyperfine ten- 
sor, the interconversion between entangled singlet and triplet 
states depends on the direction of the applied magnetic field 
through the Zeeman term in the hamiltonian. It is necessary 
also to assume that the radicals are almost immobile, without 
significant diffusive motion, in order to avoid the anisotropy 
present in the system to be averaged away (Eig. 1 (a)). The 
RP yields depend on the relative alignment of the magnetic 
field in relation to the sample [23^, 24J, so that it can work as a 
compass. 



In this work, we show that the anisotropy in the molecule 
can be replaced by an anisotropic environment (Eigure 1(b)). 
This fact allows for a free isotropic molecule in a diffusive en- 
vironment to work as a compass. In addition, although some 
discussion have been given recently about the importance of 
entanglement in the magneto-perception process 1 17, 18], this 
still remains obscure. In this sense, we also find that entangle- 
ment is neither necessary in our model of isotropic molecules 
nor in the anisotropic ones. Eurthermore, we also verify the 
functioning of the compass in the presence of artificial ra- 
dio frequency fields. We find that the isotropic model cannot 
work in the presence of such field in agreement with experi- 
mental findings |21, 25l[26l. However, the anisotropic model 
can work under some environmental conditions not unlikely 
in an open system, which notwithstanding disagree with ex- 
perimental observation. This fact gives one more evidence in 
favor of the present model of an isotropic molecule. 

Let us consider the Hamiltonian of the RP, neglecting all 
the other interactions, such as exchange and dipolar, to be: 



H = 2_^ life ' ^ik ' S/i 



h 



B S, 



(1) 



where the first term is the hyperfine contribution and the sec- 
ond one is the Zeeman contribution, with i labeling the i-th 
nucleus in the k-th radical. Earth's magnetic field is given by 
B = Bo{smO cos (/), sin sin (/), cos 0) , and A is the hyperfine 
tensor. I and S are the spin operators for the nucleus and the 
electron respectively, ge is the electron g-factor, ^o = 47/iT 
and jj^B is the Bohr Magneton. We assume that the electron 
^e -tensor is isotropic and close to that of a free electron, and 
that the hyperfine tensor A is, for simplicity, diagonal. Ad- 
ditionally, it can be either isotropic or anisotropic according 
to the model. The direction of the applied magnetic field in 











FIG. 1 . Difference between the anisotropic RPM and the isotropic 
RPM with anisotropic environment. 



relation to the RP fixed axis system is defined in terms of the 
polar angles (6>, (/)); without losing generality we are going to 
assume ^ = in order to simplify procedures. The most ap- 
propriate units to work with are those of Nuclear Magnetic 
Resonance (angular frequency units). 

In order to fully determine the dynamics of our system, we 
are going to use a master equation approach ifTvll . To account 
for the singlet or the triplet yields, we add their formation rate 
to the master equation as a dissipative process modulated by 
a factor k in angular frequency units. The simplest model we 
can consider is one nucleus coupled with its unpaired elec- 
tron plus a free electron. The evolution of the system can be 
written as 



di 



p{t)=i[p{t),H]/h 

+ ^ I {2PnP{t)Pn - PnPit) - p{t)Pn) 



(2) 



{2B^p{t)B - BB^pit) - p{t)BB^) 
{2Bp(t)B^ - B^Bp(t) - p(t)B^B). 



The operators Pn are projectors over the singlet or triplet 



shells as described in [17], and the sum is over all the differ- 
ent possible states (singlet for electronic states with nuclear 
spin up, i.e., |s, 0, t); triplet state with nuclear spin down, i.e., 
|t, 0, 1), and so on); n runs in the nuclear and electronic spin 
states. The process mediated by k can be thought as a mea- 
surement on the state of the RP, giving us information about 
the amount of singlet or triplet yield. However the density 
operator formalism enable us to take into account the envi- 
ronment in which it lives. A wide range of noise processes 
are tested to model environmental interaction which result in 
decoherence times of hundreds of microseconds ifTTll . This 
conclusion was in agreement with the experimental results 
that show an interruption of the sensitivity when a rf field 
is applied. However, it seems natural to have a faster deco- 
herence, if we take into account the very likely wildness of 
the environment surrounding the RP; we can accomplish this 
by choosing a suitable noise processes given by the opera- 
tor B. We define it as ^ = Sj^iS-2, where ^S^ |i 2} are 
the ladder operators for spins 1 and 2 respectively. The role 
of this operator is to flip spins \j3a) ^ |q^/3), leaving the 
states \aa) and |/3/3) unchanged, thus intensifying the inter- 
conversion between singlet and triplet states, and introducing 
an environment-induced asymmetry. We add this interactions 
as Lindblad-like superoperators mediated by rates P and 7. 

The yields measurement process, regardless of its magni- 
tude, is not going to affect the performance of the compass. 
Therefore k could be arbitrarily high and the compass would 
still work. However, an upper bound to k is necessary and 
it can be determined by an important experimental observa- 
tion: It was observed, in a set of experiments with European 
Robins, that an oscillating rf magnetic field Brf, perpendicu- 
lar to the Earth's one, disrupts the avian compass functioning 
lEUIJJl, leaving them without sense of direction. So when the 
magnitude of k exceeds a threshold, the influence of the mea- 
surement process should overwhelm the action of the rf field 
1 17], and as a result, the sensitivity to the variations of the 
magnetic field of Earth would be present despite of the dis- 
rupting effect. This fact can be contrasted to the experimental 
results to pick a suitable upper bound to /c as /c = O.OlMHz, 
which will be used unless stated otherwise. The processes me- 
diated by P and 7 show the same behavior, i.e., their influence 
will allow the compass to work in spite of the presence of the 
rf field if their magnitudes are high enough. As the experi- 
mental observation tell us that there will not be a compass if 
the birds are subject to this rf field, we use this to set upper 
bounds for the noise amplitudes. This fact also implies a lower 
limit for the decoherence time of the system, because an upper 
bound to the noise amplitudes implies that this time cannot be 
arbitrarily small, and we are going to have at least tens of mi- 
croseconds until the loss of all coherences in the system. We 
use a magnitude Brf = 150nT as proposed in those works. 

As stated by Schulten 1 14], in order to the compass to work 
some kind of anisotropy must be present in the system. If, 
as usual, we choose the source of anisotropy in the hyperfine 
tensor, there is a sensitivity in the RPM to all initial states. To 
show that, we start by observing that for some values of the 
rates, for example P = 7 = 2/c, the influence of the rf field 
over the compass is strengthened; and if both P and 7 are of 



^ = 7r/2 
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FIG. 2. Concurrence an Quantum Discord evolution for different 
values of the rates and magnetic field inclination angles, a) Measure- 
ment rate k = O.lMHz, and dissipation rates P = 7 = 2k, with 
and angle = 0; in the inset an angle of ^ = 7r/2 was used, b) 
Measurement rate k = O.OlMHz, and dissipation rates P = lO/c, 
7 = 2/c, with and angle ^ = 0; in the inset an angle of ^ = 7r/2 was 
used. Is worth to note that the angle affect the amount of quantum 
correlations measured by the QD and the Concurrence; for smaller 
angles the Concurrence values are not only higher, but predict the 
same time of decoherence as the QD; however, for higher angles, the 
entanglement is less important than the classical correlations, making 
the QD to be higher and to predict longer correlation times. 



the order of 20k the compass-disrupting effect by the rf field 
can still be observed. It is interesting to note that even if in 
this situation the environment contributes to the insensitivity 
of the compass by means of the rf field, in other circumstances 
it can increase the sensitivity ll27l . Also the presence of a gra- 
dient field was proven to increase the compass sensitivity f28 1 . 
One of the consequences of using the present process is that 
the decoherence times are short. To give a measure of those 
times for testing the model with anisotropy in the hyperfine 
tensor we use the Quantum Discord (QD) |29|, which is able 
to take into account a more general kind of quantum correla- 
tions, and the Concurrence |30|, which measures the amount 
of entanglement; the former is defined as 

where Iab is the mutual information between A and B and is 
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FIG. 3. Entanglement density for different values of the rates P, 
7 and /c, varying the angles of inclination of the magnetic field 
of Earth (47/iT). The entanglement was measured by the use of 
Concurrence. In the left panel, the rate values are k — O.OlMHz, 
P — ^ — lO/c, and in the right panel the values are k — 0.05MHz, 
P = lO/c, 7 = 5k. It's interesting to note that with smaller angles 
there are more entanglement sudden deaths and revivals, and that for 
angles around 7r/4 the decoherence time is shorter. 



defined as Iab = Sa-\-Sb — Sab , Sx is the von Neumann en- 
tropy of system x, J^ = max|n^} [S{pa) - Ex^^^(Pa)] 
is the classical correlation between the subsystems, Px = 
TrAiU^PABU^} and p\ = TtbIU^ PabH^J/Px', the 
maximum is taken over the positive measurements {H^} 
made over the system B. The evolution of both, Quantum Dis- 
cord and Concurrence can be seen in Figure ^. The fast loss 
of coherence is not a surprise - having an open environment 
like the one we can expect in the eye of the bird, should lead 
naturally to a fast loss of quantum correlations. This result 
should not compromise the correct behavior of the compass; 
one of the reasons, concerning the role of entanglement in the 
system, will be presented below. It is also interesting to note 
that the change in the inclination of the magnetic field also 
affects the loss of coherence. This can be seen in Figure ([3]), 
where the change of concurrence with time and inclination an- 
gle, for P and 7 fixed, is shown. For small angles there are 
going to be more collapses and revivals of entanglement, and 
for angles near 7r/4, the decoherence time are shorter that for 
angles near and 7r/2. 

Following the discussion in Ref. 1 18 1, the initial state in the 
RPM is not a perfect singlet (or triplet) state, so a natural way 
to test if the quantum correlations play a fundamental role in 
the working of the compass is to choose non-entangled ini- 
tial conditions. In fTS | several random initial conditions were 
used, and the conclusion was that there is not a crucial depen- 
dence on the entanglement of the initial state, and that even 
RPs with initial separable states with only classical correla- 
tions can produce an inclination sensitivity in the singlet yield. 
We tested several random initial conditions, some of them en- 
tangled, with the hyperfine tensor in eq. ([T]) anisotropic, and 
with the decoherence processes turned off, i.e., P = 7 = 0. 
The results show that the amount of entanglement is not a fun- 
damental factor for the sensitivity in the change of the angle 




is a must in such an open system. Both cases can be seen in 
Figure ([5]). 



FIG. 4. Angle sensitivity for non singlet or triplet initial conditions, 
specifically po = |q^<^) {oia\. Black: nuclear spin initially set to |t). 
Blue: nuclear spin initially set to a mixed state ( 1 1) + |i) )/V^- ^^^'^ 
nuclear spin initially set to ||). 



of the applied magnetic field. As an example, an initial state 
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(3) 



gives an appreciable change in the yields (and therefore allow 
sensitivity) for different angles 0. 

However, we found two cases in which without explicit 
anisotropy in the hamiltonian ([T]) there is still a fully func- 
tional compass. The first one, with P = 7 = 0, is due to 
the initial condition: if it is not an entangled one the system 
is still sensitive to changes in the inclination of the field. In 
other words, in the absence of an explicit anisotropy in the 
hyperfine tensor or in the g electronic factor, the sensitivity 
depends on the inhomogeneity of the populations in the den- 
sity matrix. As an example consider a state like eq. ([3]), i.e., 
without coherences: Its yield distribution, the blue curve in 
the Figure (|4]), shows a high variation with the angle. More- 
over, an initial state like po = |<^<^) (<^Q^|, which produces a 
lower but still appreciable angle sensitivity, generates a differ- 
ent distribution for the singlet yield depending on the nuclear 
spin state, as is shown in Figure (|4]), exemplifying the reach 
and relevance of this source of anisotropy. Now we check the 
anisotropy introduced by the environment - when there is a 
fast enough decoherence process present in the master equa- 
tion, like the processes mediated by P and 7 introduced be- 
fore. If the hamiltonian is isotropic, with singlet or triplet ini- 
tial conditions, the expected behavior is an absence of change 
in the production rates of the yields, and there is going to be 
sensitivity only if there is a decoherence process present. This 
can be understood as another class of anisotropy induced by 
the environment, which chooses a preferred direction for the 
system through the dissipation. This can open the search of a 
suitable chemical species responsible for the RP creation, be- 
cause the molecule does not need to present anisotropic hyper- 
fine or Zeeman interactions, and the degree of entanglement is 
not going to be crucial. The only requirement for the correct 
functionality of the compass is then the decoherence, which 




FIG. 5. Angle sensitivity of isotropic hamiltonians with k — 0.01. 
Black: evolution of the singlet yield for different angles with an en- 
tangled initial condition. Blue: evolution of the singlet yield for dif- 
ferent angles with an entangled initial condition and P — ^ — 2k. 
Red: Evolution of the singlet yield for different angles with an unen- 
tangled initial condition. 

Our findings go beyond establishing the role of the entan- 
glement in the compass. We found that any kind of corre- 
lation, quantum or classical, is sufficient for the RPM if the 
molecule is anisotropic, and that an isotropic molecule can 
have sensitivity for variations in the field even if it has classi- 
cally correlated, separable initial conditions. One example of 
such a state is given by eq. ([3]). Furthermore, even if they lack 
any kind of correlation but are anisotropic, giving an unbal- 
anced weight to some populations over others, like the states 
used in Figure (|4]), we can also expect a working compass. A 
careful analysis has to be made in order to transparently iden- 
tify the characteristics of the initial conditions that lead to a 
dependence of the singlet yield with the inclination angle of 
the field. 

As discussed in (231 EH, the anisotropy of a molecule 
can be averaged away if it has significant diffusive motion, 
or even rotations. Our findings of unexpected sources of 
anisotropy relax this immobility requirement. This is highly 
positive for the model given the wild conditions involving ac- 
tual bio-photochemical processes. Along with this, isotropic 
molecules seem to be more robust to environmental effects; in 
the presence of a rf field the control experimental data ll2n[25l 
shows that the compass will no longer work. However, if 
there is a strong enough dissipation the compass will work 
normally 1 17 1. An isotropic molecule can handle higher noise 
magnitudes than an anisotropic one without jeopardizing con- 
trol experimental data. 
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